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Abstract 

The mean- field theory for lossy nonlinear composites, described by complex 
and field-dependent dielectric functions, is presented. By using the spectral 
representation of linear composites with identical microstructure, we develop 
self-consistent equations for the effective response. We examine two types of 
microstructure, namely, the Maxwell-Garnett approximation and the effective 
medium approximation to illustrate the theory. 
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I. INTRODUCTION 



Recently, we proposed a mean-field theory (MFT) for nonlinear composites |IJ. The the- 
ory was applied to nonlinear transport properties of conducting composites at zero frequency, 
described by real conductivities and provided good agreement with numerical simulation [|T| . 
However, for dielectric materials near relaxation, the loss may be important. It is tempting 
to extend the MFT to complex and nonlinear permitivities. The local constitutive relation 
of the composite system is given by 

D = X (E* • E)E, (1) 

where x = x( x ) is the (position dependent) complex third-order nonlinear coefficient. In 
what follows, we denote |E| 2 = E*-E and E 2 = E-E for convenience. The effective response 
is defined as the volume average of D 

Xe E 2 E = ^ j v dVD(x), (2) 

where V is the volume of composites and E is the applied (average) field and is taken to 
be real: 

Bo = vIv dV E(X) ' (3) 
The corresponding boundary-value problem of nonlinear composite media consists of the 
field equations: 

V ■ D = 0, (4) 

V x E = 0. (5) 

There exists a potential tp so that E = V<p. The boundary condition for ip is 

if = ifo = E ■ x (6) 

on surface (S) of the composite and ipo is real. Consider the integral (1/V) J v dV D • E. 
Through the field equations, we can convert it into an integral over the surface (S) of the 
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composite: (1/V) J s dS ■ (Dip), which is equal to (1/V) J s dS ■ (Dfo) because if = <fo on 
S. Converting the surface integral back to the volume integral, we obtain: 

i f dV D • E = 1 / dV D • E . (7) 

Hence we find an equivalent expression for the effective response: 

XeK — T? f dV ~D (8) 



VJv 

Let us consider a lossy nonlinear composite in which inclusions of complex nonlinear coef- 
ficient Xii present at volume fraction pi, are randomly embedded in a host medium of \2i 
present at volume fraction p 2 , with the application of an external uniform field E . Note 
that pi + P2 = 1. From Eqs.(^) and we have therefore an expression for the effective 
response: 

XeK=PiXi(\Ei\ 2 El)+p 2X2 (\E 2 \ 2 El), (9) 

where 

(lE^E 2 ) = i ( dV |E| 2 E 2 , i=l,2. (10) 

denotes the local field average within the ith component. Similar to the derivation of Eq.([7|), 

as if* = fo on S, we arrive at 

\-\ dV D ■ E* = i / dV D ■ E , (11) 



VJv V Jv 

which implies an alternative expression for the effective response: 

Xe E* = ^((lErl 2 ) 2 ) +p 2X2 ((|E 2 | 2 ) 2 ), (12) 

where 

|2\2\ _ 1 f -7T / m?|2\2 



((|E,r) 2 ) = -/ rfV(|E| 2 ) 2 , i = l,2. (13) 
V -/Vi 

We should remark that the different types of local field averages defined in Eqs. (|TU|) and ([13]) 
are generally not equal to each other. It is interesting to note that they both give the same 



effective response [Eqs.(|9|) and (|12|)1 . However, the nonlinear partial differential equations 
pertaining to the boundary-value problem cannot be solved analytically. We shall invoke 
the mean-field theory [[L]] to obtain an approximate expression for the effective response. 
According to previous work JTJ, we consider the linear constitutive relation: 

D = eE, (14) 

where e = e(x) is the (position dependent) complex dielectric constant. By definition, we 
find the effective linear response: 

e e E = — f dV eE. (15) 



VJv 

From Eq. ([?]), we have therefore 

e e El = p 1 e 1 (Ej}+p 2 e 2 (El}. (16) 

From Eq.(|ll|), we have alternatively 

e e E 2 = Pl e 1 (\E 1 \ 2 )+p 2 e 2 (\E 2 \ 2 ). (17) 

This expression will be useful for calculating (|Ej| 2 ). In order to estimate Xe, we invoke the 
decoupling approximation 

(|E 4 | 2 E, 2 ) = (|E J | 2 )(E 4 2 ), i = l,2. (18) 

within the ith component. This assumption is good in microstructure for which the local 
electric field is nearly uniform within the ith component, but less accurate when these 
fluctuations are large, as in a random composite near the percolation threshold ||. With 
this assumption, 

Xe E^ = PlXl (|E 1 | 2 )(E 2 )+p 2X2 (|E 2 | 2 )(E 2 ). (19) 

Comparing Eqs. (|T6"D and (JT9]), it is tempting to write ei = XidE^ 2 ), e 2 = x 2 (|E 2 | 2 ) and 
e e = XeEg. Hence, within the decoupling approximation, we may interpret the nonlinear 



component as a linear dielectric material with a field-dependent dielectric constant. Suppose 
e e is known as a function of its constituent dielectric constants, 

e e = F(e 1 ,e 2 ,p 1 ), (20) 

then we find the local field averages |jj 

Pl (El) = ^K (21) 
P2 (El) = ||Eg. (22) 

However, we have to determine (|Ej| 2 ) in the ith component, which is generally not equal 
to (Eif). Thus a straightforward application of the previous formalism |E] is impossible. We 
resort to an alternative approach based on the spectral representation 0. 

In what follows, we shall consider two important types of microstructures: (1) Dispersion 
microstructures as in the Maxwell-Garnett approximation [[| and (2) symmetric microstruc- 
tures as in the Bruggeman effective medium approximation ||. 

II. MAXWELL-GARNETT APPROXIMATION 

The Maxwell-Garnett approximation (MGA) is good for dispersion microstructures and 
the theory is inherently non-symmetrical ||. For convenience, we consider the case in which 
component 1 is embedded in component 2. For a linear composite of E\ and e 2 at volume 
fractions p\ and p 2 respectively, the MGA reads ||]: 

e e — e 2 ei — e 2 

Pi ~H — T\ ' ( 23 ) 



e e + (d-l)e 2 r ei + (d-l)e 2 

where d is the dimensionality of composites. Solving Eq.(|23|), we obtain 

_ E 2 [ei(dpi + p 2 ) +e 2 (d- l)p 2 ] 
P2E1 + (d -p 2 )e 2 

From Eqs.(pl|) and (p2|), we can calculate the local field averages: 



(24) 



Pid 2 e\ 
[p 2 ei + (d- p 2 )e 2 ] 



pi<e?> = — r)v\,, 2 K (25) 
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,21 



P 2 (E ) = 1 - ^\ J E . (26) 

In MGA, the local field in the inclusion is uniform, we can determine Ex and hence E^ 
explicitly, we find 

Pi (|Ei | 2 ) = P f^\ , 2 Eg. (27) 

\p 2 ei + [d-p2)e 2 \ 2 

From Eq.(|I7|) and using Eqs.(^) and (|27j), we find 



P 2 (|E 2 | 2 ) = 1 - \ ,s \ ' J E o- (28) 



We should remark that although |(E 2 )| = (|Ei| 2 ) in MGA, |(E|)| is generally not equal 
to (|E 2 | 2 ). In fact, |(E|)| < (|E 2 | 2 ). A common error in the literature is to treat these 
averages equal even when x is complex. If we write e\ = xi(\^>i\ 2 ) and e 2 = X2(|E 2 | 2 ), then 



Pid[d\e 2 


2 -£> 2 |ei - e 2 


2 ] 


|P2£i - 


f (d- 


P2)e 2 
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Eqs.(|27j)-(28) can be solved self-consistently for (|Ei| ) and (|E 2 | ) and hence the effective 
nonlinear response can be calculated. 

In Fig.|I|, we present the MGA results in three dimensions (3D). We let xi = 1 + 3« 
and X2 = 3 + i and compute Xe as a function of p\. As pi increases, Re(x e ) decreases from 
Refe) = 3 towards Re(xi) = 1 while Im(x e ) increases from Im(% 2 ) = 1 towards lm(xi) = 3. 
Since the contrast between the two components is relatively small, the local field averages 
( | Ei | 2 ) and (|E 2 | 2 ) remain close to unity. However, unlike |(E 2 )|, |(E|)I shows a significant 
deviation from (|E 2 | 2 ), indicating the nonuniformity of local field in the host. 

For two-component composites, it has proved convenient to adopt the spectral represen- 
tation of the effective linear response ||: Let v = 1 — e 1 /e 2 , w = 1 — e e /e 2 , and s = 1/v, we 
find 

, . f 1 m(s')ds' 

where m(s) is the spectral density which is obtained through a limiting process: 

m(s) — lim Im w (s + irf) . (30) 

Furthermore, m(s) obeys the sum rule 



/ m(s')ds' = p\. (31) 
Jo 

Eqs.(|24])-(|2~8"|) can readily be converted into the spectral representation. We find 

w = (32) 
s — Si 

where Si = P2/d. From Eg. (|30|) , we obtain the spectral density in MGA: 

m(s) = p\5(s — si). (33) 
Similarly, the local field averages [Eqs. (|25|)— (p8|)l are given by 

< e ?> - 4?' (34) 



^ = jCT < 35 > 



M^=(l- P -^=^£)K (36) 



(1-V3i) 

P2 (|E 2 R=(l-^igil)E S . (37) 

If we let v = 1 — xi(|Ei| 2 )/x2<|E 2 | 2 > and solve Eqs.(gD and © self-consistently for <|E! | 2 > 
and (IE2I 2 ), the effective nonlinear response can be determined. Although the solution 
appears somewhat simpler, the same results are obtained. 

III. EFFECTIVE MEDIUM APPROXIMATION 

The Bruggeman effective-medium approximation (EMA) is known to be symmetrical 
with respect to interchanging components 1 and 2 ||. Note that there is a percolation 
threshold in the theory. For a random linear composite of ei and e 2 at volume fractions p± 
and P2 respectively, the self-consistency equations reads 0: 

ei — e e e 2 — e e . . 

Pi ~ , , — 7\ + P2 - j-, — tt— = 0, 38 
ei + (d-l)e e e 2 + {d-l)e e 

where d is the dimensionality of composites. Solving Eq. (|3"8"D , we obtain e e and hence w and 

the spectral density 



m{s) = ^-^8{s)B{d Pl - 1) + mi (s), (39) 

where S(s) and 9(d P i — 1) denote the Dirac delta function and the Heavyside step function 
respectively, and 
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m i( s ) = ^ 71 — TVv( s ~ Sl )( s 2 - s), s 1 <s<s 2 , (40) 
Z7i[a — L)s 

= 0, otherwise, 
where Si and s 2 are given by 

si = \ (l + (d - 2) Pl - 2^(d-l)p 1 (l-p 1 )) , 



s 2 = - M + (d - 2) Pl + 2^{d - l)p x {\ - Pl )j . (41) 

When pi > 1/d, which is the percolation threshold in EMA, component 1 percolates the 
composite and there is a 5 function contribution to m(s) at s = 0. Eqs.(|3"4])-(|3"TD can readily 
be generalized to ||: 

,2v f 1 m(s'W ^ 2 



"^ = 1 jr^ E - < 42) 



Pl <| El [ 2 > = / pr^^Eg, (43) 



n/F 2^ A f 1 (l-v 2 s')m(s')ds' \ 2 

P2 (E 2 ) = ^1 - ^ _ — j E , (44) 

p 2<IE2 | 2 >^ri-r (i -':' 2 ^r )d ^ E g- w 



One can check that when the MGA spectral density [Eq.(|33"D1 is used, Eqs.(|34])-(|37l) are 
recovered. Again we let v — 1 — % x { |Ex | 2 )/ X2 ( |E 2 1 2 ) - Eqs. (|4"3|) and ([4"5|) are coupled integral 
equations. They can be solved numerically for (|Ei| 2 ) and ( | E 2 1 2 ) and hence the effective 
nonlinear response can be determined. 

In Fig.[| we present the EMA results in 3D. We let xi = 1 + 3« and \2 = 3 + i. FigfJ 
exhibits similar behavior as Fig.l, however, with some differences. Here Re(xe) and Im(x e ) 
cross at pi = 0.5 because EMA describes the symmetric microstructure and we have the 
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symmetric choice of xi an d Xi- Other quantities, like the local field averages (|Ei| 2 ), (|E 2 | 2 ), 
and |(E 2 )|, |(E 2 )| also cross at pi = 0.5. In EMA, both the local fields in the inclusion and 
host regions are not uniform. 

DISCUSSION AND CONCLUSION 

The model may be applicable to lossy dielectric composite materials under an intense 
electric, available from laser. The present formalism can naturally be extended to optical 
properties of nonlinear composites [[/J. When the size of the inclusion is smaller than the 
wavelength of light, we can use the quasi-static limit. The mean-field theory can be used to 
calculate the optical nonlinearity. 

In conclusion, we extend the mean-field theory for lossy nonlinear composites, described 
by complex and field-dependent dielectric functions. We developed self-consistent equations 
with the aid of the spectral representation of the linear composites. We examine two impor- 
tant types of microstructure, namely, the Maxwell-Garnett approximation and the effective 
medium approximation to illustrate the present theory. 
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FIGURES 

FIG. 1. Real and imaginary parts of the effective nonlinear response plotted against volume 
fraction of the inclusion in the MGA with %i = 1 + 3i and X2 = 3 + i. 

FIG. 2. Real and imaginary parts of the effective nonlinear response plotted against volume 
fraction of the inclusion in the EMA with %i = 1 + 3z and \2 = 3 + i. 



11 



MGA 




o L ^— — L ^ ' ^ ' ^ ' ^ 1 

0.0 0.2 0.4 0.6 0.8 1.0 

P 

EMA 




1 1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 



